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The fluid theory of plasmas is extended to include the properties of electron spin. The 
linear theory of waves in a magnetized plasma is presented, and it is shown that the spin 
effects causes a change of the magnetic permeability. Furthemore, by changing the direction 
of the external magnetic field, the magnetic permability may become negative. This leads 
to instabilities in the long wavelength regimes. If these can be controlled, however, the spin 
plasma becomes a metamaterial for a broad range of frequencies, i.e. above the ion cyclotron 
frequency but below the electron cyclotron frequency. The consequences of our results are 
discussed. 

I. INTRODUCTION 



The theoretical concept of materials with a negative index of refraction has a long history, while 
the practical developments are relatively recent in this field. In his Moscow school, Mandelstam 
presented results_for materials exhibiting a negative index of refraction (see also the historical 



account in Ref. 



3; 



Negat™ pernnttMty and permeabil.ty materials, and the eonsecuenees of 

a negative group velocity, were moreover given attention by Pafomov [4] and Agranovich and 
Ginzburg 5|, and have since then been discussed by several authors, most notably by Veselagojsl], 
who gave a detailed analysis of the consequences of such material properties (see also 0, 1?,!^, 
Although not known to be found naturally, such materials have recently been realized in laboratory 
environments [K 



(see Refs. 



12 



and the experimental development in conjunction with theoretical insi ghts 



13|) has spawned a rapidly growing interest in these materials (see, e.g. 
for a review). From a sharp resonance in the material response to the applied external field, one 
may obtain negative e and /i. The normal procedure for obtaining negative- index-of-refraction 
materials is to put together two structured materials that both have negative permittivity and 
neg ative permeability, such that the resulting composite material has a negative refractive index 



16l | . A nonlinear metamaterial can also be constructed through the nonlinear properties of the 



constituent materials 



13, 



18t | , admitting new types of solitary wave structures 



19|, 



2C 



2li 



The field of quantum plasmas is a rapidly growing field of research. From the non-relativistic 
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domain, with its basic description in terms of the Schrodinger equation, to the strongly relativistic 
regime, with its natural connection to quantum field theory, quantum plasma physics provides 
promises of highly interesting and important application, fundamental connections between differ- 
ent areas of science, as well as difficult challenges from a computational perspective. The necessity 
to thoroughly understand such plasmas motivates a reductive principle of research, for which we 
successively build more complex models based on previous results. The simplest lower order effect 
due to relativistic quantum mechanics is the introduction of spin, and as such thus provides a first 
step towards a partial description of relativistic quantum plasmas. 

Already in the 1960's, Pines studied the excitation spectrum of quantum plasmas 



22, 



231], for 



which we have a high density and a low temperature as compared to normal 

interest in the 



there has 



31 



been an increased 



32 



33, 



34 



35 



36 



371, 



38 



3S 



40 



aroperties of quantum plasmas 



development in nanostructured materials 
plasmas 



42, 



431, 



24 



plasmas. Rece ntly 



25 



26 



27|, 



4J]. The studies has been motivated by the 



471, 



45( 1 and quantum wells 



461 ] ■ the discovery of ultracold 



48|, or a general theoretical interest. The list of quantum mechanical effects that can 



be included in a fluid picture includes the dispersive 



)otential 



24 



25|, 



effects 



24 



26 
49 



25 



27 
50 



26 



27 



28 



29 
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28 
51 



spin properties 



31 



37|, 



32 



381, 



33 



39|, 



34 



id, 



aarticle properties accounted for by the Bohm 



41 



52] • Within such descriptions. 



35 



361 ]. the zero temperature Fermi pressure 



24 



as well as certain quantum electrodynamical 



251, 



261, 



27 



37, 



38|, 



5d, 



51 



521 ] quantum 



and classical collective effects can be described within a unified picture, sometimes even showing a 



4l|. 



surprising overlap between classical and quantum behaviour 

Here we study the linear theory of electromagnetic wave propagation in a magnetized plasma, 
with a special focus on the properties caused by the electron spin. We are then able to present a 
scheme for such a system to display metamaterial behavior. Specifically this is induced by exposing 
a low temperature high density plasma to an external magnetic field, which creates a magnetization 
in the plasma due to the electron spin. By changing the direction of the external magnetic field, the 
magnetic permeability may become negative. It should be noted that the above procedure induces 
instabilities in the long wavelength regime. A number of ways to control these instabilities are 
pointed out. Assuming that this can be done successfully, the spin plasma becomes a metamaterial 
for a broad range of frequencies, i.e. above the ion cyclotron frequency but below the electron 
cyclotron frequency. The conditions needed to create a sufficient magnetization is discussed in the 
final section of the manuscript. 
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II. BASIC EQUATIONS 



The theory for 



291, 
spin 



3d 



31 



37 



32 



38|, 



33 



quantum 



34 



35 



plasmas including the effects o " 



36l |. the Fermi pressure 



24 



25 



26 



partic 



271, 



e dispersion 



24j, 



25 



26 



27 



281] and effects due to the electron 



39[ has been described in a number of recent papers. For our purposes it will be 



sufficient to include the spin effects, as the Fermi pressure and the Bohm-De Broigle potential will 
not affect whether the plasma is a metamaterial or not.^ Furthermore, neglecting terms that are 
quadratic in the spin vector,^ the governing spin plasma equations can be written 



dt 



+ V • (n,v,) = 0, 



(1) 



rusus [ol + ^s- Vj V, = (E + V, X B) + -Tls^VBa - VP,, (2) 

together with the spin evolution equation 

| + v,.v)s=5^1s.xB, (3) 

Here Qs and is the charge and mass of species s, /u, = qsh/2ms is the magnetic moment, Ps = 
ksTVug is the pressure(for simplicity we use an isothermal pressure model), is the number 
density, is the velocity and is the spin vector. 

Next we concentrate on the linear wave modes in a magnetized plasma described by Eqs. ([ID-(I3D 
together with Maxwell's equations 

VxE = -- (4) 

and 

I dE 

VxB=^oJ+^^ (5) 
where, in addition to the free current density, we include the magnetization current source 

= V X M, = V X (^^s,) , (6) 



For this purpose we chose a coordinate system such that the unperturbed magnetic field is Bg = 
Bqz, the wavevector is k = k±x + /c^z, and the variables are divided into an equilibrium value 



^ This is clear from the dispersionrelation presented in [52], where the effect of the Fermi pressure and the Bohm de 
Broigle potential is incorpatated in an effective thermal velocity, as far as the linear wave modes are concerned. 

^ The terms that are quadratic in the spin vector can be omitted if the characteristic spatial scale is longer than the 
thermal de-Broigle wavelength. 
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(index 0) and a perturbed part (index 1). For simplicity index 1 are omitted on the electric field 
and the velocity, since these variables have a zero unperturbed part. As a preparation, we first 
consider the linear theory without the spin terms. Using the continuity equation to express the 
density in terms of the velocity, the momentum equation relates the velocity to the electric field. 

Then from the momentum equation we solve for the velocity in terms of the electric field to find 
the susceptibility tensor for each particle species. Combining this with Maxwell's equations, the 
result becomes 



Sij ( 1 



ijj'- 



+ 



+ Xii = 



(7) 



with bij is the Kronecker delta and the susceptibility tensor is 



-klvl 



\ 



i 5 



,u)?,,u}csk±kzv?. 



uij^k±kzvf^ 



u)^^u)csk±kzv^^ 



(8) 



where vts = {ksT /rngf'^'^ is the thermal velocity, Ups = {noQ^ / £oitt-s)^^^ the plasma frequency, 
ojcs = qsBo/ms the gyrofrequency and uj\ = (w^ — k^v"^^ (a;^ — k'^v'^^) ^^cs {^"^ ~ ^z^L) "^l^z'^ts- 
The next aim is to add the electron spin contribution, where the spin effects due to the ions is 
neglected due to their small magnetic moment. In general the spin vector is a dynamical variable 
whose relation to the EM-field is complex already in the linear theory. Thus from now on, we limit 
ourselves to the case where the dynamics is slow compared to the spin precession period, which is 
equivalent to limiting ourselves to frequencies lo <C \uJce\ In that case we can write the (electron) 
spin vector as 



s = 4Btanhf^ 
2 V ksT 



(9) 



where the tanh-factor is due to thermodynamic considerations.^ The linearized magnetization 
current then becomes 



= tanh ^Vni x z + noV 



Bn 



(10) 



^ Using Fermi-dirac statistics, the slight overweight of particles with spin orientation in the lower energy state results 
in a macroscopic spin- vector proportional to tanh. (nBBo/kBT). 
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The magnetic field Bi can be expressed in terms of E from Faradays law, and the density is 
given in terms of the velocity from the continuity equation, which is expressed in terms of E from 
the given susceptibility of each species. Thus the magnetization current (|10p is also expressed in 
terms of E. However, before that procedure is implemented, we must also modify the standard 
susceptibility Eq. ([8]) to account for the spin dependent force in the momentum equation. This can 
be achieved by noting that when solving the momentum equation for the velocity, the magnetic 
dipole spin force can simply be accounted for by including the different components as "effective 
electric fields" . Thus when solving for the electron velocity in terms of the electric fields, the spin 
force is included simply by making the substitutions 

E,. = £;. + tanh(^) ^Bi, 

E,y = Ey (11) 

Again expressing Biz in terms of E through Faraday's law, these alterations can be expressed as 
a spin modification of the free current susceptibility. Thus formally we can write 

/ = ffree + 4 = X^ee^j + X%Ej 

where the direct spin magnetization contained in xip can be determined from (jlOp (by expressing 
Bi and ni in terms of E), and the free part of the susceptibility can be divided as 

i'j ij I ij 

Xfree Xl ' Xmd 

where is the (free current) susceptibility due to the Lorentz force given by dSD, and the contri- 
bution from the magnetic dipole force x^d ^^''^ be found from ^ combined with the substitution in 
(jlip and the z-component of Q. The theory outlined here is straightforward, but results in rather 
cumbersome formulas. To reduce the complexity, and arrive at more transparent expressions, we 
introduce the following simplifications: 

1. The plasma is quasi- neutral, which is a valid approximation provided uJ^^ S> 

2. The perpendicular (to the magnetic field) free electron (i.e. non-spin) part of the current can 
be approximated by the E x B-drift, which is valid for frequencies well below the electron 
gyro frequency. 

3. The displacement current in Maxwells equations is small. This is valid when point 1 applies 
together with uj'i^ » a;^, and amounts to neglecting the Kronecker delta term in ([7]) 
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4. Only electron thermal motion is of significance, which is valid if k'^v}- <C uj'^. 



The theory outlined above now reduces to Eq. ([7]) with Xij given by 



\ 



Z^i 



'Ti 



.i^i.ik\k,ci 



(12) 



where ujpi = ujpi{l — ^qMq/ Bq)^^/"^ is the spin modified ion plasma frequency, = [c^ — 
[^BBo/mi)ta.VLh.{^BBQ/kBT)Y/'^ is the spin modified ion-acoustic velocity and = vti{me/miY/'^ 
is the standard ion-acoustic velocity. Furthermore, the unperturbed magnetization Mq is given by 
Mq = riQfiB tanh.{fiBBo/kBT). We note that by contrast to the other tensor components, X33 is not 
modified by the spin effects. As a specific example, we can investigate the magnetohydrodynamic 
limit uj <^ u>ci, in which case we get the dispersion relation 



LO' - klCl 



UJ 



k'Cl 



kll»A[u:^-klcl)-klk 



,2 ,,2~4 



0. 



(13) 



by putting the determinant of ([7]) to zero. Here Ca = Ca{^ — fJ-oMo/Bo)^/'^ is the spin-modified 
Alfven velocity, and Ca = (Bq/ fiominQ)^^"^ is the standard Alfven velocity. Eq. (fT3]) agrees with 
the results of Refs. \4Ul\ in the appropriate limiting cases, provided a sign error is corrected in the 
last term of their dispersion relation. The first factor of Eq. (|13|) describes the shear Alfven mode, 
whereas the second factor has two roots, describing the fast and slow magnetosonic modes. 



III. THE POSSIBILITY OF A PLASMA AS A METAMATERIAL 



Next we are interested in the possibilities to get metamaterial properties. As a starting point 
we note from (fT2|) and ([7]) that the (relative) dielectric component £xx is given by 

,2 



(jJ. 



pi 



(14) 



which become negative for frequencies above the ion cycclotron frequency. Here the first Kronecker 
delta term in ([7]) has been neglected due to point 3 above. Next we note that the relative magnetic 
permeability is given by /i,. = Bq/{Bq — /ig-Mo), which causes the transition from w^- — > uj^^ 
in (jl2p . In thermodynamic equilibrium, the spin magnetization Mq = ^u^no tanh(/iBi?o/A;Br) 
enhances the external field, which correspond to a paramagnetic situation where //^ > 1- However, 
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let us now assume that we have a laboratory plasma immersed in an external magnetic field, where 
the internal spin magnetization gives a significant contribution to Bq. Then consider what happens 
if we rapidly switch direction of the external field 180°, and study the properties before the spin 
state has time to reach a new thermodynamic equilibrium state. In particular we are interested 
in the case where the external contribution to Bq is smaller than the internal contribution ^jlqMq, 
and directed in the opposite direction. The above linearized theory then still applies, but with 
the difference that fir = Bq/{Bq — ^qMq) is now is a negative quantity. A first observation of the 
changed properties of this system is that long wavelength waves (i.e. kCA < ^d) described in (jl3p 
is now unstable, since C\ changes sign with (jl^ and becomes negative for the above scenario. We 
shall assume that it is still useful to study the stable waves with shorter wavelengths, however. In 
particular this is of interest in case one of the following conditions apply 

1. The plasma system is of a rather small size, and the long wavelength waves ( i.e. UCa < 
ojci) are stabilized by inhomogeneities not included in the model. 

2. The growth rate 7 of the long wavelength (at most of the order 7 ~ i^ci) is slow enough such 
that there is still time to study the physics on a time scale much shorter than a;~^ 

3. Some dissipative mechanism not included in our model will be sufficient to stabilize the 
instability (however, note that the addition of a finite resistivity will not suffice for this 
purpose). 

As a specific example we study transverse waves with to > Ud propagating along Bq (that is in 
the opposite direction of Bq — ^qMq). The dispersion relation then reduces to: 



k'Ci (15) 



where the +(— ) sign corresponds to a right hand (left hand) circular polarization. As described 
above, both fir and changes sign with {Bq — ^qMq), and thus ([15]) confirms that waves are 



Wa 



> 2u) cii the waves 



unstable in the long wavelength regime. However, for shorter wavelengths, 
are stable independently of the sign of C\. Next considering this special case, such that k and uj 
are real, we find that the time-averaged Poynting vector (S) is 

(S) = (ExH) = |E|^^ ^'^°"^°^°^ (16) 

UJ Bq 

which apparently changes sign with Bq— ^qMq such that the above system shows the characteristics 
of a metamaterial. 
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IV. SUMMARY AND CONCLUSION 



Previous theories 



371, 



38( 1 on magnetized spin plasmas have been extended. In particular, the 
linear theory have been generalized to cover the frequency range tOci ~ w ^ ^ce- The purpose 
have been to investigate whether it is possible to produce the characteristics of a metamaterial. It 
is found that in principle this is possible by switching the direction of an external magnetic field 
180°, provided the spin magnetization of the plasma is sufficiently large. However, it is a great 
challenge to produce the desired plasma conditions in the laboratory. In particular we need to 
combine low temperatures with high densities to obtain a sufficient magnetization. Ways to reach 
high plasma densities have been known for a rather long time, and methods to reach extremely low 



plasma temperatures has recently been found 



47 



48( 1 ■ On the other hand, the ultra cold plasmas 



are still of too low density to be useful for the above purpose. Laser produced plasmas can reach 
sufficient densities, and it is possible that experiments can be designed to keep the temperature 
sufficiently low. A challenge with such a setting might be that the plasma background dynamics 
is too fast for successful experiments of this kind to be done. In conclusion, the production of 
a plasma metamaterial consists of two challenges. Firstly, to produce a plasma with a sufficient 
magnetization, and secondly to master the long wavelength instabilities that are introduced when 
the direction of the external field is changed. 
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